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ABSTRACT 

We present a semi-analytic, physically motivated model for dark matter halo concen¬ 
tration as a function of halo mass and redshift. The semi-analytic model combines an 
analytic model for the halo mass accretion history (MAH), based on extended Press 
Schechter (EPS) theory, with an empirical relation between concentration and forma¬ 
tion time obtained through fits to the results of numerical simulations. Because the 
semi-analytic model is based on EPS theory, it can be applied to wide ranges in mass, 
redshift and cosmology. The resulting concentration-mass (c — M) relations are found 
to agree with the simulations, and because the model applies only to relaxed halos, 
they do not exhibit the upturn at high masses or high redshifts found by some recent 
works. We predict a change of slope in the z ~ 0 c — M relation at a mass scale of 
10 11 Mq. We find that this is due to the change in the functional form of the halo 
MAH, which goes from being dominated by an exponential (for high-mass halos) to 
a power-law (for low-mass halos). During the latter phase, the core radius remains 
approximately constant, and the concentration grows due to the drop of the back¬ 
ground density. We also analyse how the c — M relation predicted by this work affects 
the power produced by dark matter annihilation, finding that at z = 0 the power 
is two orders of magnitude lower than that obtained from extrapolating best-fitting 
c — M relations. We provide fits to the c — M relations as well as numerical routines 
to compute concentrations and MAHB 
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1 INTRODUCTION 

Over the past few years large cosmological simulations have 
been performed to determine the properties of dark matter 
halos, includ ing density profiles, shapes and accretion histo¬ 
ries ( see e.g. ISpringel|l2005l : iKlypin et al.ll201ll : iBrvan et alj 
l2013h . These properties are of particular interest, as forming 
galaxies depend on the structural properties of the halos in 
which they are embedded. 

During hierarchical growth, halos acquire a density pro¬ 
file with a near universal shape, that can be described b y a 
simp le formula known as the ‘NFW profile’ (I Navarro et alj 
1 199 A hereafter NFW). The NFW density profile is described 
by just two parameters, halo mass, M, and concentration, 
c. A halo’s concentration is defined as the ratio of the virial 
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radius, R v n, and the scale radius, r_ 2 , which is defined as 
the radius where the logarithmic density slope is —2. Thus, 
given the NFW profile, only a relation between concentra¬ 
tion and halo mass (hereafter, the c — M relation) is needed 
to fully specify halo structure at fixed mass. Therefore, nu¬ 
merous studies have been undertaken to improve the c — M 
calibration. 


Despite its importance, there is still no solid agree¬ 
ment on the dependence of halo concentration on halo mass 
and redshift. A small change in the adopted cosmology can 
hav e important effects on the structure of dark matter ha¬ 
los (jMacdo_et_alJ 2008]), and on their mass accretion histo¬ 
ries ( Zhao et aljl200Sf) . For example, the mean concentra¬ 
tions of dwarf-scale dark matter halos change by a factor 
of 1.5 between the various W ilkinson Microwave Aniso tropy 
Probe (WMAP) co smologies dSper gel et al .l2003l.lg007l) . The 
Planck cosmology dPlanck Collaboration fc et alj 2014] ') has 
higher matter density, H m , and higher power spectrum nor¬ 
malization, erg, compared to the cosmologic al parameters of 
the year 5 data release of WMAP fWMAP5: lKomatsu et al] 
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l2009f) . The Planck cosmology therefore suggests that halos 
assemble e arlier and are more conc entr ated (c.f. c— M re la- 
tions from lDutton fc Maccidll2014l and lDuffv et al.l|2008h . 

However, cosmology may not be the primary reason for 
the differences in the concentration—mass relations found by 
various authors. Recent works that adopt the same cosmol¬ 
ogy still_^nd_differentc— M relations (compare for exa m¬ 


ple Pntfxm_fc_Maccib 2014 and Picmcr & Kravtsov 20151. or 


IKly pin et al.l 201 ll and Prada et al.ll2012h . Putton fe Maccibl 


( 2014ll found that the c—M relation is well described by a 


power law, but flattens at high r edshift and exhibits a posi- 
tive slope at 2 > 4. In contrast, IPiemer fe Kravtsovl (l2015il 
found a strong upturn in the high-mass end of the c—M 
relation at all redshifts. The disparity between these stud¬ 
ies could be due to the dy namical state of the selected dark 
matter halos. For example. iLudlow et al.l J2012ll showed that 
massive halos that are substantially out of equilibrium are 
more likely to be found at a transient stage of high concen¬ 
tration, thus explaining the puzzling upturn in the high-mass 
end of the c—M relation. Indeed, they reported that the 
upturn disappears when onl y dynamic ally-relaxed systems 


are considered. However, [Kl^|3in_et_cT (|2014) argued that 
by“ 


the virial criterion used bv lLudlow et alj ( 2012 ) to select re¬ 


laxed systems is incorrect, as it needs to include effects o f 
the surface pressure and external forces. iKlypin et al.l (|2014l ) 
modified the virial criterion and ended up selecting massive 
halos t hat had previ o usly b een considered as unrelaxed. As a 
result. IKlypin et alj (1201414 obtained an upturn in the c—M 
relation of their relaxed halo sample and concluded that the 
upturn is a real feature of the c—M relation. They explained 
that as extremely massive halos have more radial infall ve¬ 
locities, infalling mass penetrates deeper within the inner 
halo, thus increasing the concentration and producing the 
upturn. 

The main goal of this work is to derive a physically 
motivated model for the c—M relation of relaxed halos based 
on the dark matter halo accretion history. By relating the 
concentration to the halo accretion history, we find that the 
c—M relation does not show any upturn or strong flattening. 
We then study the c — M relation in detail using simulations 
and selecting relaxed halos without using the virial criterion, 
and investigate whether recently accreted particles are able 
to reach the inner parts of the halo and thus increase the 
concentration. 

Our c — M model relies on the fact that concentrations 
depend on the evolutionary stage of halos when they were 
formed. Several works have suggested that halo formation 
can be described as an ‘inside out’ process, where a bound 
core (of a certain fraction of the halo mass today) collapses, 


logical accretion rate (iManriaue et al. 20031: 

Wang & White 

20091: Dalai et al. 2010: Salvador-Sole et al. 

201211. In this 


framework, the halo concentration should depend on the 
epoch at which a certain fraction of the halo mass was as- 


sembled. As a result , various authors Bullock_et_aLl 1 200 
IWechsler et al.|[2002l : IZhao et al.l 120031 : Ludlow et ali bi 
have provided mo dels that relate c to the halo mass his¬ 
tory. For instance, IZhao et al.l (|2003h showed that when the 
mass accretion rate of a halo slows down at low redshift, 
its scale radius, r_ 2 , remains approximately constant, and 
hence that concentration scales with the virial radius. On 
the other hand, in the regime of a high mass accretion rate 


(at high redshift), the scale radius scales approximately as 
the virial radius and thus c remains constant. 

The connection between a halo’s mass accretion his¬ 
tory (hereafter MAH) and its concentration, c, is there¬ 
fore obtained through its ‘formation’ time. The halo forma¬ 
tion (or assembly) time is traditionally defined as the point 
in time when the halo mass reached a fraction of the to¬ 
tal mass today. Low-mass halos typically assemble earlier, 
when the Universe was denser, than high-mass halos do. As 
a result, low-mass halos are more concentrated. Clearly, if 
concentration correlates with formation time, and forma¬ 
tion time depends on the mass variance, a (because a de- 
scrib es the halo MAH, se e the analytic model for the MAH 
from lCorrea et al.ll2015al ). then it is expected that c corre¬ 
lates with a and hence with the peak height, v, defined as 
v = 1.686/cr. This is indeed what several works have found 

fe.g. IZh ao et al.| [2j )Qgl:[Prada et al1l2012l:lL udlow et al.||201l|: 

iDutton fe Maccibl 2014 ). We showed in Correa et al. ( 2015bf) 

that the physical origin of the c — a (or c—v ) relation is the 
halo MAH. 

Re cently, in ICorrea et al] (l2015all and ICorrea et al.l 
d2015bh (hereafter Paper I and Paper II, respectively), we 
provided two models for the MAH of halos, an analytic 
model and a semi-analytic model. The semi-analytic model 
uses a functional form for the MAH, that is motivated by 
extended Press-Schechter (EPS) theory, and links the MAH 
to halo structure through two empirical relations obtained 
from simulations. The analytic model is fully derived from 
the EPS formalism and thus does not require calibration 
against any simulation data. 

In the analytic model, the halo MAH is described 
in terms of the rms of the density perturbation field, 
er, as M(z) = Mo(l + where Mo refers 

to the present halo mass, a depends on cosmology, and 
/(Mo) ~ 1/cr(Mo). This expression illustrates that as a de¬ 
creases with halo mass, the function /(Mo) increases, caus¬ 
ing the exponential in M(z) to dominate. As a result, high- 
mass halos accrete faster than low-mass ones, due to their 
low value of a. As low a implies large peak height, the EPS 
formalism predicts that density perturbations with large v 
experience an accelerated collapse phase relative to the av¬ 
erage, and grow faster in time. 

In this work we present a semi-analytic, physically mo¬ 
tivated model for dark matter halo concentration as a func¬ 
tion of halo mass, redshift and cosmol og y. T h e sem i-analytic 
model, which builds on that of lLudlow et alj (12014I) . uses the 
analytic model for the halo MAH provided in Paper I, as well 
as an extension of the empirical relation between concentra¬ 
tion and formation time obtained through fits to simulations 
provided in Paper II. As a result, the semi-analytic model 
for halo concentrations shows how the c—M relation is ex¬ 
pected to evolve based on the hierarchical growth of halos. 

This paper is organized as follows. We begin in Section[2] 
with a description of the set of cosmological simulations used 
in this work. In Section [3] we describe the analytic MAH 
model provided in Paper I and extend it to high redshift. 
In Section [5] we define halo formation time and build an 
empirical relation between formation time and concentra¬ 
tion through fits to simulation data. Next, we describe the 
semi-analytic model for halo concentrations that combines 
the analytic model for the MAH and the empirical relation 
described previously. We analyse the evolution of concentra- 
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tion that predicts the semi-analytic model in Section [5] In 
Section [G] we discuss the impact of the results of our semi- 
analytic model for halo concentration on the signal from 
dark matter annihilation. In Section [3 we discuss the main 
assumptions the semi-analytic model relies on. Finally, we 
summarize and conclude in Section [8] 


2 SIMULATIONS 

Throughout this work we compare our analytic results to 
the output from numerical simulations. We use a set of cos¬ 
mological dark mat ter o nly (DMON LY) simulations from 
the OWLS project (jSchave et all 120101 ). These simulations 
were run with a significantly extended version of th e N-Body 
Tree- PM, SPH code gadget3 (last described in ISprineell 
12005t h The initial conditions were generated with CMB- 
FAST (version 4.1; ISeliak fe Zaldarriagalll996l ~) and evolved 
to redshif t 2 = 127, wher e the simulations were started, 
using the IZerdoviclJ fl97flh approximation from an initial 
glass-like state Jwhitelll99al . In order to assess the numer¬ 
ical convergence, we use simulations of different box sizes 
(ranging from 25 to 400 /i -1 Mpc) and particle numbers 
(ranging from 128 3 to 512 3 ). The simulation names con¬ 
tain strings of the form LxxxNyyy, where xxx is the sim¬ 
ulation box size in comoving h _1 Mpc and yyy is the cube 
root of the number of particles. Our DMONLY simulations 
assume the WMAP5 cosmology. However, to investigate the 
dependence on the adopted cosmology, we use an extra set 
of five dark matter only simulations (100h -1 Mpc box size 
and 512 3 dark matter particles) which assume values for the 
cosmological parameters derived from different releases of 
the WMAP and the Planck missions. See the tables in Ap- 
pendix[X]for the sets of cosmological parameters adopted in 
the different simulations, as well as the main numerical pa¬ 
rameters of the runs such as comoving box size, number of 
dark matter particles, dark matter particle mass, comoving 
gravitational softening and maximum physical softening. 


3 HALO MASS ACCRETION HISTORY 

We begin this section by briefly reviewing the analytic model 
for the MAH derived from the EPS formalism in Paper I, 
and showing how the MAH depends on cosmology and on 
the initial peak of the primordial density field. In Section 
we extend it to estimate the halo MAH tracked from 
an arbitrary redshift. Readers only interested in the c — M 
relation model can skip directly to Section [4] 

3.1 Analytic model for the halo mass history 

In Paper I, we used simple analytic arguments based 
on the E PS forma lism and the analytic formulation of 
iNeistein et alj J2006l l. to show that the ‘shape’ of the MAH is 
determined by the growth factor of the initial density pertu- 
bations. The halo MAH is well described by an exponential 
in the high-redshift regime, but it slows to a power law at 
low redshift, because the growth of density perturbations is 
halted in the dark energy dominated era due to the acceler¬ 
ated expansion of the Universe. Therefore, we showed that 
the expression 
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Figure 1. Halo MAH of a 10 12 M© halo (coloured lines) obtained 
from the model given by eqs. and by assuming various 

cosmologies as indicated in the legend. The grey lines correspond 
to MAH obtained from DMONLY simulations that assume the 
Planck and WMAP5 cosmologies. In the top right corner, we plot 
a versus halo mass, to show that the change in a under different 
cosmologies drives the change in the MAH. 


M(z) acdm = Mo(l + z) a e^“, (1) 

accurately captures the median halo MAH, where Mo refers 
to halo mass today, and a and /3 are parameters that de¬ 
pend on Mo, cosmology and the linear power spectrum. In 
the case of an Einstein de Sitter (EdS) cosmology (S7a = 0 
and = 1) or an open universe (Da = 0 and < 
1), there is no acceleration in the expansion of the Uni¬ 
verse at low redshift. Then the halo mass history is sim¬ 
ply described by an exponential as M{z )gas = Moe^ z , where 
/3 = -~1.686(2/7r) 1,/2 /(Mo). For a complete description of the 
model, see Paper I. 

We find that the MAH model can be used to calculate 
halo mass histories in cosmologies other than WMAP5, and 
that the differences are mainly driven by the changes in as 
and fl m . We show this in Fig. [l] where the halo MAH of a 
10 12 Mg halo (coloured lines) was estimated for the various 
cosmologies, as indicated in the legend. In the top right cor¬ 
ner of Fig. [T] we plot a versus halo mass, to show how the 
change in a drives the change in the MAH. The exception 
is the Planck cosmology, which has a relatively low as but a 
large = 0.317, which raises M(z) close to the WMAP1 
M(z). 

The overplotted grey lines in Fig. [Q correspond to the 
MAH obtained from DMONLY simulations that assume the 
Planck and WMAP5 cosmologies. In this case, we compute 
the MAH of the main subhalo (that is not embe dded inside a 
large r halo) of Friends-of-Friends (FoF) groups dDavis et al.l 
Il985ll . by tracking the virial mass of the main progenitor at 
each prior output redshift. Halo virial masses and radii were 
determined using a sp herical overdensity ro utine within the 
SUBFIND algorithm (ISpringel et al.ll200lh centred on the 
main subhalo of FoF halos. Throughout this work we define 
the halo mass as the total mass within the radius r 2 oo for 
which the mean internal density is 200 times the critical 
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Figure 2. Median MAHs for halos of 10 11 Mg starting from various redshifts. In both panels the grey solid lines correspond to the 
analytic model described in Section 13.21 The coloured curves in the left panel are the MAHs obtained from the DMONLY simulations 
WMAP5—L025N512 and WMAP5—L050N512. The mass histories are computed by calculating the median value and the la error bars 
are determined by bootstrap resampling the hal os from the merger tree at a given output redshift. The coloured dot dashed curves in 
the right panel are the MAHs obtained from the Ivan den Bosch et alJ +014 ) model. 


density. For a more detailed description of the method used 
to create merger trees, see Paper II. 

3.2 Analytic model for the MAH: high redshift 
prediction 

The model presented in Paper I is suitable for estimating 
halo MAHs that are tracked from z = 0. In this section we 
extend this analytic model to estimate MAHs of halos of the 
same halo mass that are tracked from arbitrary redshifts Zi. 
This is shown in Fig. [2] where the MAHs of 10 11 Mq halos 
are obtained from DMONLY simulations (coloured curves). 
The curves show the mean MAH of halos of the same mass 
(10 11 Mq in this case) that begin at Zi = 0 (blue curve), 1 
(dark green curve), 2 (green curve), 3 (orange curve) and 
4 (red curve). High-redshift MAHs are dominated by large 
accretion rates and characterized by a pure exponential. 

We generalize the analytic model so that it describes the 
MAHs from any z% redshift. Expression JlJ can be rewritten 
as 

M(z,M(zi),Zi) = M(zi)( 1 + z — Zi) a e^ z ~ zi \ (2) 

where M(z, Zi) denotes the MAH of a halo with mass M(zi) 
at redshift Zi. In the above expression, z > Zi and the pa¬ 
rameters a and /3 depend on M(zi) and redshift Zi 


a = 


( 3 ) 

P = 

-f(M( Zi )), 

( 4 ) 

f(M(zi)) = 

[a 2 (M{zi) / q) - a 2 (M(zi)))~ 1/2 , 

( 5 ) 

a\R) = 

^ P(k)W 2 (k; R)k 2 dk, 

(6) 

<1 = 

4.137 x z f -°' 9476 , 

( 7 ) 

Zi = 

—0.0064(lo glo Mo) 2 + 0.0237(lo glo M 0 ) 
+1.8837, 

(8) 


where D(z) is the linear growth factor, P(k) the linear power 
spectrum, W (fc; R) the Fourier transform of a top hat win¬ 
dow function and R defines a in a sphere of mass M = 
(47r/3)p m ,of? 3 , where p mj o is the m ean background density 
today. We use the approximation of lEisenstein Hul d 19981 1 
to compute P(k), normalized such that a(8/i _1 Mpc) = as- 
As a result, /(Mo) depends on the power spectrum and 
halo mass. It can be seen from eqs. © and © that 
at large Zi, a —> 0 due to D(zi) oc 1.686(2/7r) 1 ^ 2 /(l + 
Zi) for Zi 1, indicating that the MAH is mainly de¬ 
scribed by an exponential. Table [I] provides a summary 
of the nomenclature adopted throughout this work. The 
above equations introduce an analytic halo MAH model 
directly derived from EPS theory that does not require 
calibration against any simulation data (see Paper I for 
more details). The numerical values given in eqs. © and 
© were determined by assuming the WMAP5 cosmology 
(n m = 0.258, Ha = 0.742, h = 0.72, n s = 0.963, ag = 0.796). 

In the left panel of Fig. 0 we compare the model given 
by eqs. ©-© to various MAHs obtained from a set of 
DMONLY simulations. Our analytic model is shown by 
grey solid lines, where we have taken M(zi ) = 10 11 Mq. The 
coloured curves in the left panel correspond to the MAHs ob¬ 
tained from the DMONLY simulations WMAP5-L025N512 
and WMAP5-L050N512. We find very good agreement be¬ 
tween the simulation outputs and the analytic model at 
all redshifts. The simulation outputs from the boxes L = 
25h _1 Mpc and L = 50/i _1 Mpc converge up to 3 = 5. At 
higher z, the outputs from the L = 25h~ 1 Mpc simulation 
underestimate M(z) because the box size limits the maxi¬ 
mum sizes of the structures that can form at each redshift. 

In the right panel of Fig . [2) we compare our e xtended 
analytic model with t he Ivan den Bosch et al.l (l2014l l model. 
Ivan den Bosch et alJ ll2014l) e xtracted halo mass histories 
from the Bolshoi simulations (iKlypin et alJ 1201 il l and ex¬ 
tended them below the numerical resolution limit using EPS 
merger trees. Once they had obtained the MAH curves for 
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a large range of redshifts and halo masses, they made use of 
a semi-analytic model to transform the (average or median) 
MAHs, based on the Bolshoi cosmology, to other cosmolo¬ 
gies. Using their publicly available code, we calculate the 
mass history curves for the WMAP5 cosmology for compar¬ 
ison with our results. We find that ther e is some discrepancy 
at hig h-redshift for all the curves. The Ivan den Bosch et all 
(12014 1 MAH model seems to over predict the halo mass at 
z > 5, most likely as a consequence of the different halo def¬ 
initions, and subtle differences in the definition of the main 
progenitor (van den Bosch, private communication). Overall, 
there is very good agreement between the most recent accre¬ 
tion history study in the literature and our model, as well 
as with the simulation outputs. In Section [4] we will make 
use of our analytic MAH model to calculate concentrations. 

Using the extended MAH model for high redshift, we 
can calculate the accretion rate of a halo at redshift z. We 
differentiate eq. © with respect to time and replace dz/dt 
by — Ho[n m (l + z ) 5 + Ha(1 + z) 2 ] 1 ^ 2 , to obtain 


Table 1. Notation reference. 


Notation 

Definition 

M[z{) 

Total halo mass at Zj, 
defined as M 200 

r 200 or r 2 00 [M(zi),Zi] 

Virial radius at Zi of a halo of 
total mass M(zi) 

r -2 or r_2[M(zj), Zi\ 

NFW scale radius at Zi 

c or c[M(zi), Zi\ 

NFW concentration at Zi 

M r (r , Zi) 

M(< r), mass enclosed within r 
at Zi of a halo of total mass M(zi) 

M r (r_2,Zi) 

Mass enclosed within r_ 2 at Z{ 

M(z, M{zi), Zi) 

Mass at 2 : of a halo with mass 

M(z{) at Zi 

Z-2 

Formation redshift, when equating 
M(z- 2 ,M(zi),Zi) to M r (r- 2 ,Zi) 

(p)(< r-2,Zi) 

Mean density within r _2 at Zi 

Pcrit,0 

Critical density today 

Pcrit { z i) 

Critical density at Zi 


dM(z, M(zi), Zi) —i (M{z,M(zi),Zi)\ 

- - - = 71.6M 0 yr ^ 1 O ! 2 M0 ) 

x (^ 7 ) [-a/(l + z-zi)-j3] 

x (1 + z) [$l m (l + z) 3 + SIa] 1 / 2 , (9) 

where a and $ are given by eqs. © and 0, respectively. 
Note that the above formula will give the accretion rate at 
redshift z of a halo that has mass M(zi) at redshift z», and 
mass M(z, M(zi), Zi) at redshift z. 


4 CONCENTRATION- MASS RELATION 


A theoretical understanding of the physical connection 
between concentration (the parameter that characterizes 
the internal structure of NFW dark matter halos) and 
the initial conditions of the density field, is essential 
for the physical interpretation of relations like c — u 
(concentration—peak height) or c — M, that have been cal- 


200ll: iNeto et a.1.1120071: Gao et al 

120051: Maccio et all 20071: 

Duffv et al. 200sl: Ludlow et al. 

20131: Dutton & Maccio 

20141: Diemer & Kravtsov 2015|). 


It has previously been shown that concentration is de¬ 
termined by the halo MAH, and that the MAH depends on 
the p ower spectrum and the adopted cosmo l ogical param- 
eters jWechsler et al.l 120021 : IZhao et al.l 120031 : ILudlow et al.l 
120131 . l2014lf . In this section we show, through analytic and 
numerical modelling, how the concentration of dark matter 
halos depends on cosmology and the power spectrum of den¬ 
sity perturbations. Our results imply that the halo MAH is 
the physical link between concentration and peak height. 


4.1 Formation redshift 


halo’s MAH and its internal structure is the formation red¬ 
shift. For a halo with mass M(z;) at redshift Zi, we define 
the formation redshift to be z_ 2 , the redshift at which the 
mass of the main progenitor equals the mass enclosed within 
the scale radius at z = z», 

Z —2 = z[M(z-2, M(zi) , Zi) = Af r (r_ 2 ,z*)], (10) 

dLudlow et al.ll2013h . Here M(z_ 2 ,M(zj), Zi) is the mass at 
Z— 2 of a halo with mass M(zt) at Zi, and we denote the mass 
enclosed within r, M(< r), as M r . For an NFW profile the 
internal mass M r {r- 2 ,z%) is related to the total halo mass 
as 


M r (r- 2 , Zi) 


M(%) 


Y( 1) 

Y(c[M(zi),Zi]) ’ 


( 11 ) 


where Y(u) = ln(l+u)— u/(l+u), c[M(zi),Zi] is the concen¬ 
tration at Zi and M(zi) is the total halo mass at Zi. In cases 
where we identify halos at Zi = 0 and track their mass histo¬ 
ries, we calculate z_2 by setting M(z~ 2) equal to the mass 
enclosed within r_2 today. In cases where we identify halos 
at Zi > 0, we first calculate r_2 and M r (r_2,Zi) at the par¬ 
ticular redshift Zi, and then find Z-2 by tracking the MAH 
(for z > Zi) and equating M(z_2, M(z;), Zi) to M r (r- 2 , Zi). 
See Table [J for a summary of the nomenclature. 

ILudlow et al.1 (l2013l . l2014l) and Paper II showed that z_2 
correlates strongly with c, and in Paper II we demonstrated 
that the scatter in z_2 and in the halo MAH predicts the 
scatter in c. In this section we explore how the formation 
time — concentration relation varies for halos identified at 
various redshifts. 

We computed density profiles and MAHs for halos iden¬ 
tified at redshifts z. ( = 0,1,2,3 and 4. The density profiles 
were computed by fitting the NFW density profile, 


0 ( r ,.\ = _ PcTit(zj)5 c _ 

PK, * } (cr/r 2 oo)[l + (cr/r 2 oo)] 2 ’ 


( 12 ) 


As discussed in the Introduction, halo MAHs can be used 
to estimate halo concentrations. Halo concentrations re¬ 
flect the mea n density of the Universe at the formation 
redshift (NFWlWechsler et al.ll2003 : IZhao et aT1l2003l . 120091 : 
ILudlow et al .11201 .'lli . Therefore, the essential link between a 


for each individual halo. In the above equation, 
Pcrit(-Zi) = 3H 2 (zi)/8nG is the critical density of the 
universe, 5 C is a dimensionless parameter related to the 
concentration c = r2oo/r~2 by 5 C = and >"200 is the 

virial radius. 
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Figure 3. Left panel: mean density within the NFW scaled radius, (p)(< r_ 2 ), at Z{ against the critical density of the universe at 
the formation time, p C rit(-z— 2 )- Each dot in the panel corresponds to an individual relaxed halo identified at z% and coloured by mass 
according to the colour bar at the top of the plot. The open symbols show the median values of the sample at z% as indicated by the 
legend in logarithmic mass bins of width 51og 10 M = 0.4. The solid line shows the best linear fit to the p c rit(^— 2 ) — (p)(< 2 ) relation. 

Right panel: formation redshift against concentration. The solid lines show the c — z —2 relations given by equation GU for various z-i as 
indicated in the legend of the left panel. The open symbols correspond to the median values of the samples in logarithmic mass bins of 
width 5 log |q M = 0.4 and are colour coded by z r . The grey areas show the scatter in z_ 2 . 
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Figure 4. Concentration—mass relations at z = 0 (left panel), z = 1 (middle panel) and z = 2 (right panel). The dots in the 
panels correspond to individual, relaxed halos illustrating the scatter in the relation. The simulations assume the WMAP5 cosmological 
parameters and have box sizes of 400, 200, 100, 50 and 25 h~ x Mpc, as indicated. Because of resolution limits only halos in the mass 
ranges indicated in the top legend were used from a particular simulation. The open symbols show the median c—M values in logarithmic 
mass bins of width <5 log 10 M = 0.4. The solid line shows the prediction of the c — M model obtained from the halo MAH as described 
in Section [I] 


We begin by fitting NFW profiles to all halos at Zi that 
contain at least 10 4 dark matter particles within the virial 
radius. Throughout this work, we define the virial radius as 
7 *200, the radius for which the mean internal density is 200 
times the critical density. Then, for each halo, all particles in 
the range —1.25 ^ log 10 (r/r 2 oo) 0 are binned radially in 
equally spaced logarithmic bins of size A log 10 r = —0.078. 
The density profile is then fitted to these bins by perform¬ 
ing a least-squares minimization of the difference between 


the logarithmic densities of the model and the data, assum¬ 
ing equal weighting. The corresponding mean enclosed mass, 
M r (r-2,Zi), and mean inner density at r_ 2 , (p)(< r_ 2 ,z»), 
are found by interpolating along the cumulative mass and 
density profiles from r = 0 to r _2 = r 200 /c, where c is the 
concentration from the fit of the NFW halo. Then we gener¬ 
ate merger trees for these halos and by interpolation we de¬ 
termine the redshift Z-2 at which M{z- 2 , Zi) = M r (r- 2 , Zi). 













In order to obtain robust estimates and to test 
whether the c — M relation includes an upturn in th e me¬ 
dian concentrations of massive halos (IPrada et al.l 120121 : 
iDutton fe MacciDl2014l : iDiemer fe Kravtsovll2015l iT we only 
consider ‘relaxed’ halos. We define relaxed halos as those 
halos for which the separation between the most bound par¬ 
ticle and the centr e of mass of the FoF halos is smaller t han 
0.07-Rvi r (fo llowing iMaccio et al.ll2007l . iNeto et al.ll2007l and 
iDuffv et al.ll2008h , where is the radius within whi ch the 
mean density is A, as given bv lBrvan fe Normanlll998l . times 
the critical density. Our relaxed sample contains 2425 halos 
at z = 0, 726 halos at z = 1, 226 halos at z = 2 and 78 and 
20 halos at z = 3 and z = 4, respectively. 

The left panel of Fig. [3] shows the mean density within 
the NFW scale radius, r_2, at redshift Zi. The median values 
of (p}(< r_2,Zi) follow the best-fitting relation 


(p){< r - 2 , 


Zi) = 200 


c[M{ Zi ),ZifY( 1) 
Y(c[M(zi), Zi]) 


Pcrit ( Zi ), 


(13) 


expressed as a function of the critical density of the Universe 
at Z— 2 , 


Pcrit ( 2 — 2 ) — Pcrit,o[fIm(l + Z— 2 ) T I2a] , (14) 

where p cr it,o = 3 H 2 (z = 0)/8nG. Note that densities along 
both the x— and y —axes are expressed in units of the crit¬ 
ical density at Zi. Each dot in the panel corresponds to an 
individual relaxed halo identified at z; and coloured by mass 
according to the colour bar at the top of the plot. The open 
symbols show the median value of the sample in logarithmic 
mass bins of width <flog 1(1 M = 0.4 and are coloured by Zi 
as indicated in the legenqj. 

At each redshift Zi, the p cr i t (z_ 2 ) — (p)(< r_ 2 , Zi) corre¬ 
lation clearly shows that halos which collapsed earlier have 
denser cores. 

We perform a least-squares minimization of the quantity 
A 2 = j, Ef=iKPt>(< r- 2 ,Zi) - i ? (p crit ,j(z_ 2 ), A)], where j 
goes from 1 to the number of dark matter halos, N, at Zi and 
F’(p cr it,j(z_ 2 ), A) = A x Pcrit,j(z— 2 ), to obtain the constant 
of proportionality, A. The solid line corresponds to the best- 
fit to the Pcrit (z — 2 ) — (p)(< r_ 2 , Zi) relation, and we find (in 
agreement with Ludlow et al. 2014) that the average relation 


(P)(< r- 2 ,Zj) Pcrit (*•'— 2 ) , . 

Pcrit (.Zi) Pcrit (.Zi) 1 ' 

is maintained through time with A = 887 ± 36, where the 
la error was obtained from the least-squares fit. 

Using eqs. m and GD we can rewrite this relation as 


c[M(zi),zi] 3 y(l) A [flm(l T Z—Yy’ T Ha] , . 

Y (c[M(zi), Zi]) 200 [On(l +Zi) 3 + 0\] ’ 1 ’ 

The right panel of Fig. [3| shows the c — z ~2 relation (solid 
lines) given by eq. (fl6|) for various Zi. The open symbols 
correspond to the median values of the sample in logarithmic 
mass bins of width 5 log 10 M = 0.4. The grey areas show the 
scatter in z- 2 . 


1 Note that it is possible for individual halos to appear multiple 
times in Fig. \3 \ (left panel). For example, a 10 13 Mq halo at z = 0, 
has a total mass of ~ 10 12 ' 2 Mq at z = 2, therefore the halo will 
be included in the p cr it( z — 2 ) — (p)(< r — 2 ,Zi) relation at Zi = 0 
but also at Zi = 2. 
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4.2 Semi-analytic model for halo concentration 


In this section we describe the semi-analytic model for halo 
concentration as a function of halo mass and redshift. This 
model combines the analytic model for the halo MAH given 
by eqs. m and the empirical relation between z _2 and c 
given by eq. m- 

We begin by calculating M{z- 2, M(zi), Zi) from eq. ([2]), 
and use the equality 


M(z- 2, M(zi), Zi) M r (r-2, Zi) T(l) 

M(zi) M(zi) Y(c[M(zi),Zi])’ ( 

which follows from eqs. m and m and is valid under the 
assumption that the halo density profile follows the NFW 
profile, to obtain 


Y(l) 

Y(c[M(zi),Zi]) 


(1 + Z-2 - 


Zife^ z - 2 ~ zi) , 


(18) 


where a and /3 are given by eqs. m and Q, respectively. 
Next, we combine eqs. Gi and m to obtain the concentra¬ 
tion, c[M(zi), Zi], of a halo of total mass M(zi) at z». We re¬ 
mind the reader that throughout this work the adopted halo 
mass definition is M 200 , and the concentrations are therefore 
defined as c = C 200 ■ 

Fig. [4] shows the concentration—mass relation at z = 0 
(left panel), at z = 1 (middle panel), and at z = 2 (right 
panel). The dots in the panels correspond to individual re¬ 
laxed halos identified in the simulations at Zi = 0, 1 and 
2, whereas the open symbols correspond to the median val¬ 
ues in logarithmic mass bins of width <51og 10 M = 0.4. The 
solid line shows the c—M relation that results from the semi- 
analytic model described above. We find excellent agreement 
between the median values from the simulations and the 
c — M relation predicted by the semi-analytic model at all 
redshifts. 

So far we have adopted the WMAP5 cosmology. In Ap¬ 
pendix [B] we discuss the dependence of our c — M relation 
model on cosmology and extend it to make it suitable for 
any values of the cosmological parameters. 


4.3 Impact of relaxedness on the c—M relation 

Several recent studies llK lvpin et al.ll201ll; Pradaetal ] 2oIi 
IDutton fe Maccidl 120141 : Diemer fe Kravtsov 2015l i have 
found that the c — M relation flattens at high redshift and 
exhibits an ‘upturn’ at the high-mass end, meaning that the 
concentration increases with halo mass for the most massive 
halos. In this section we investigate whether this interest¬ 
ing behaviour is seen in our semi-analytic model or in the 
simulation outputs. 

Our model does not predict an upturn. The model re¬ 
lates c to the MAH via the formation redshift, z_ 2 (see 
Fig. [3] right panel), which decreases with halo mass, mean¬ 
ing that more massive halos are less concentrated because 
they formed more recently. If c were to increase with halo 
mass, then high-mass halos would have to form earlier than 
low-mass ones, at a point when the Universe was denser. 
This behaviour is neither seen in our simulations (see Fig.[T] 
coloured lines), as we only consider relaxed systems, nor pre¬ 
dicted by EPS theory, because it would be antihierarchical 
for such systems. 

To investigate further, we use the simulation outputs 
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Relaxed halo sample Complete halo sample 





Figure 5. Concentration — mass relation at 2 = 0,1, 2, 3 and 4 under the WMAP5 cosmology for the relaxed halo sample (left panel) and 
the full sample (middle panel). The open symbols indicate the median concentrations in logarithmic mass bins of width <51og 10 M = 0.25 
at z = 0 and z = 1, and <5 log 10 M = 0.30, 0.40 and 0.50 at z = 2, 3 and 4, respectively. Only bins containing at least 10 halos are shown. 
The error bars show the la confidence limits. The dashed lines correspond to the best-fitting power laws to the open symbols. In the left 
panel, the solid grey line shows the c — M relation predicted by the semi-analytic model. The right panel shows the fraction of relaxed 
halos, with respect to the complete sample, for each mass bin and redshift. The inclusion of unrelaxed halos results in a flattening of, or 
even an upturn in, the c — M relation at high redshift. 
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Figur e 6. Comparison of the c — M relation predicte d by our model (solid lin es) with the relations of van den Bosch et al.| [ 20 1 4) (1 eft- 
panell. fPiemer fc Kravtsov] il20 1 Oil ) (middle-panel) and lPutton fc Macciol il20 1 4l) (right-panel). Note that PuttorT' fo'^accicJ 120 1 1) consider 
only r elax ed halos for their best-fittin g model and fit their model in the halo mass range 10 10 — 10 15 M/.,. whereas |plemer"fe Kravtsov! 
( 2015 1 and lvan den Bosch et a! ] (203) use their full halo sample and fit their model in halo mass range 10 10 —10 15 M@ and 10 11 —10 15 Mg, 
respectively. 


to calculate concentrations by fitting NFW profiles to lralos 
that are resolved with at least 10 4 particles wit hin the viria l 
radiu s, and for which the convergence radiufl (jPower et al.l 
120031) is smaller than the minimum fit radius of 0.05 times 
the virial radius. In addition, we consider two halo samples. 


2 The convergence radius is defined such that the two—body dy¬ 
namical relaxation time-scale of the par ticles is simi l ar to the 
age of the univers e. For more details see I Power et alJ l2003ll or 
IPuffv et al.l ll2008l l. 


A relaxed halo sampltjfl and a full halo sample. When con¬ 
sidering only relaxed halos, as we have done so far, we find 
that we restrict our halo sample to around 80% of the total 
at 2 = 0, 65% at 2 = 1, 55% at 2 = 2, 50% at 2 = 3 and 
43% at 2 = 4. 

Fig. [5] shows the c — M relation (at various z) of the 


3 As proposed by iNeto et al.l (12007! ) , relaxed halos are defined 
as those halos for which the separation between the most bound 
particle and the centre of mass of the FoF halo is smaller than 
0.07 times the virial radius. 
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relaxed sample (left panel) and the full sample (middle 
panel). These panels show the median value of the concen¬ 
tration (open symbols) in logarithmic mass bins of width 
51og 10 M = 0.25 at 2 = 0 and 1, and 51og 10 M = 0.30, 0.40 
and 0.50 at 2 = 2, 3 and 4, respectively. We increase the bin 
size with redshift so that each bin at a fixed mass contains on 
average approximately the same number of halos. For each 
bin the lcr error bars were determined by bootstrap resam¬ 
pling the halos. Only bins containing at least 10 halos are 
shown. The dashed lines correspond to the best-fitting power 
laws to the open symbols. In addition, the left panel shows 
the c — M relations predicted by the semi-analytic model in 
the solid grey lines. The middle panel shows a strong flatten¬ 
ing a nd upturn in the c— M relati on at high 2 , in agreemen t 
with iMuhoz-Cuartas et ahl J201ll ') and IPrada et all (|2012l h 
However, this upturn is not seen for the relaxed sample. 
Thus, we conclude that the previously seen upturn results 
from the inclusion of unrelaxed halos, in agreement with 
iLudlow et al.l J2012l 'l. We show the fraction of relaxed halos 
(with respect to the total sample) for each mass bin and red- 
shift in the right panel of Fig. [o] We find that the relaxed 
fraction tends to decrease towards high mass and redshift. 

Our results suggest that the dynamical state of dark 
matter halos should be considered when analysing the pa¬ 
rameters that describe the halo internal structure, because 
the density profiles of unrelaxed halos are poorly captured 
by the NFW fitting formula fe.g. iNeto et allf2007h . Because 
halo concentrations are clearly affected by transient depar¬ 
tures from equilibrium, we only consider relaxed halos in the 
remainder of this work. 


4.4 Comparison with previous studies 

In this section we compare the c — M relations of 
the most re cent studies on da rk matt er halo con- 
va^_de^^osch_et_alJ_ (l2014l l (hereafter, 


centrations, 


vdB14), iDiem er fe Kravtsov! ( 2015li (hereafter, DK14), 
Dutton fe Macciol ( 2014ll (hereafter, DM14) and 


Ludlow et a l. ( 2014! ) (hereafter, L14), with the model 

presented in this work. _ 

vdB14 used the c—M relation of lZhao et al.l (|2009l l (ob¬ 
tained from fits of a full halo sample from numerical simu¬ 
lations) and adjusted the parameters by fitting it to the 
c — M relation of the full halo sample from the Bolshoi sim¬ 
ulations. vdB1 4 assumed the Bolshoi c osmology (consistent 
with WMAP7. iKomatsu fe et al.ll201 il l, but they made use 
of a semi-analytic model to scale their model to any cosmol¬ 
ogy. We assume the Planck cosmology and use the publicly 
available code of vdB14 to calculate their c—M relation. 
DK14 obtained a concentration model given by a best-fit 
seven parameter function of peak height ( v ) and slope of 
the linear power spectrum. They considered their full halo 
sample and extended their model to make it suitable for any 
cosmology. Finally, DM14 followed the evolution of the con¬ 
centration of relaxed dark matter halos from a series of N- 
body simulations that assumed the Planck cosmology. DM14 
fitted a power-law to the c — M relation and restricted their 
analysis to relaxed halos only. 

The left panel of Fig.[ 6 ]shows a comparison of our c — M 
model (solid lines) to the model of vdB14. To compare with 
vdB14, we predicted the concentrations using the analytic 
expression for the MAH assuming the Planck cosmology 


(shown in Fig. [lj and a 2-2 — c relation with a constant 
of proportionality of 850 instead of the value 887 used for 
the WMAP5 cosmology (see Appendix [B] for a discussion 
of the cosmology dependence of our model). We find broad 
agreement with the relation of vd B14 only at 2 — 1 and 2. In 
their w ork. Ivan den Bosch et al.l (l2014l l used the lZhao et al.l 
d2009l l model which assumes that c never drops below 4 at 
high redshift. 

The middle-panel of Fig.[ 6 ]shows the DK14 c — M rela¬ 
tion calculated assuming the Planck cosmology. As they in¬ 
cluded their entire sample of halos for their c — v relations, 
they obtained an upturn at the high-mass end at all red- 
shifts. We find that our model predicts concentrations that 
are a factor of 1.2 larger just before the high-mass upturn. 
Finally, the right-panel of Fig. [ 6 ] shows reasonable agree¬ 
ment between our model and the DM14 c—M relation for 
z = 0,1,2 and 3 although the results diverge at low masses. 
In their work, DM14 fitted a power-law, c oc M a , to the 
c — M relation at all redshifts, and found that the slope, 
a, increases from —0.1 at z = 0, to 0.03 at 2 = 5. As they 
restricted their halo sample to relaxed halos, they did not 
obtain a significant upturn at the high-mass end of the c—M 
relation. 

Fig. [ 6 ] shows that the physically motivated model pre¬ 
sented in this work yields c — M relations that are generally 
in agreement with previous results. However, the important 
improvement with respect to previous works is that we are 
presenting a physical analytic model that can then be ex¬ 
trapolated to very low-masses, and is suitable for any cos¬ 
mology. 

The model for dark matter halo concentrations pre¬ 
sented in this work strongly relies on the relation (p(r_ 2 )) — 
Pcrit(2— 2 ), which supports the idea th at halos grow inside- 
out. This relation was introduced in ILudlow et al.l (120131 1 
and explored in L14, who recently presented a related 
model for the concentration-mass r elation. In thei r wor k, 
L14 used the a v erage MAHs from Ivan den Boschl (120021 4 
and IZhao et al.l d2009h that begin at Zi = (Q. They fit¬ 
ted the halo MAHs, written as M(p cr it), with the NFW 
profile expressed in terms of the enclosed density. They 
looked for a correlation between the concentration parame¬ 
ter cmah, that results from an NFW fit to the halo MAH, 
and the concentration parameter of the halos density pro¬ 
file, cnfw, and used the best-fitting relation to predict 
cnfw from cmah- L14 and this work use the same forma¬ 
tion redshift definition to connect concentrations with halo 
MAHs. L14 used the (p(r_ 2 )) — p c rit( 2 - 2 ) relation to find 
the cmah — cnfw relation, whereas in this work we used 
the analytic MAH model to define formation redshifts and 
used the (p(r_ 2 )) — p cr it(.z- 2 ) relation to predict concentra¬ 
tions. Although there is good agreement between L14 and 
our c—M relation at 2 = 0, there are differences in the 
relations at high— 2 , e.g. a factor of 1.2 difference between 
the concentrations of a lO lo /i - 1 M 0 halo at z = 2 (c ~ 5.25 
versus cli4 ~ 6.3), and a factor of 1.58 for a lO s /i _ 1 M 0 halo 
at 2 = 2 (c ~ 7.95 versus cli 4 ~ 12.58), for the WMAP5 


4 L14 c — M mod el used MAHs from Ivan den Boscbl <|2002| 4 and 
IZhao et al.1 DoOfl ) to show specific examples on how to construct 
a c(M, z) relation for a given MAH, but any MAH model can be 
used. 
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Figure 7. Predicted concentration — mass relation for the WMAP5 cosmology over a wide range of halo masses (log | 0 Ml Mq = [—2,16]) 
and redshifts (z = 0 — 20). The solid lines correspond to our c — M model obtained from the halo MAH as described in Section [I] The 
lines are coloured as a function of redshift as indicated. The dashed lines correspond to power-law fits to the high-mass c — M relation. 


cosmology. Those differences are mainly due to the different 
MAH models. Since the (p(r_ 2 )) — p C rit( 2 - 2 ) relation is es¬ 
sentially equivalent to the cmah — cnfw relation, we expect 
L14 and our semi-analytic model to give consistent results 
if the same MAH model is used. We believe however that 
we have improved upon the L14 c — M model by combin¬ 
ing the (p(r_ 2 , zi)) — p C rit(z— 2 , Zi) relation with an analytical 
MAH model, M(z,Zi), that begins at any redshift Zi, and 
allows a detailed analysis of the redshift dependence of the 
c(M, z) relation for relaxed halos. Another important differ¬ 
ence is the tentative evidence for a cosmology dependence 
in the (p(r-2,Zi)) — p a it(z-2, Zi) relation (for a discussion 
see Appendix ®. 


4.5 Extrapolation to low halo masses and high 
redshifts 

Because our semi-analytic model for halo concentration is 
physical, rather than a purely empirical fit to the simula¬ 
tion results, we can use it to extrapolate beyond the mass 
and redshift ranges spanned by our simulations, assuming 
that the z_ 2 — c relation given by eq. (fl6l) holds. Fig. [T] 
shows the predicted concentration-mass relation for a wide 
range of halo masses (log 10 M/Mg = [—2,16]) and redshifts 
(z = 0 — 20). The dashed lines correspond to the high-mass 
power-law c—M relations at low redshift. These are included 


to aid the comparison of the slopes of the c — M relation in 
the high- and low-mass regimes. There is a clear ‘break’ in 
the z = 0 c — M relation. For M > 10 12 Mq concentration 
scales as c oc m~°' 033 , whereas at M < 10 9 Mq it scales as 
c oc M~°' 036 . The change of slope around these halo masses 
is substantial up to z = 3 — 4. However, at z > 4 there is 
no ‘break’ in the c — M relation. In Section [5] we provide a 
tentative explanation for the physical origin of the break in 
the c — M relation. 

We provide fitting functions for the c—M relation in 
the high-z and low-z regimes. The following expression is 
suitable for the low-redshift regime (z ^ 4) and at all halo 
masses, 

logi 0 c = a + 0 log 10 (M/ M 0 )[l + 7 (log 10 M/ M 0 ) 2 ], (19) 
a = 1.62774 - 0.2458(1 + z) + 0.01716(1 +z) 2 , 

/3 = 1.66079 + 0.00359(1 + z) - 1.6901(1 + z ) 0 00417 , 

7 = -0.02049 + 0.0253(1 +z ) -0 1044 . 

In the high-redshift regime the c — M relation can be 
fitted using only two parameters. The following expression 
is suitable for z > 4 and at all halo masses, 

logi 0 c = a + P log 10 (M/ M 0 ), (20) 

a = 1.226 -0.1009(1 + z) +0.00378(1 + z) 2 , 
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Figure 8. Top panel : evolution of the concentrations of halos that 
at 2 : = 0 have masses of Mq = 10 6 ,10 8 ,10 10 ,10 12 and 10 14 Mq, 
as indicated in the legends. Bottom panel : MAHs of halos of the 
same masses as in the top panel. 


/3 = 0.008634 - 0.08814(1 + 2) -0 ' 58816 


The above fitting functions have been calculated assum¬ 
ing the WMAP5 cosmology. Appendix[B]provides a series of 
best-fitting relations for the Planck cosmology, as well as a 
short discussion of the cosmology dependence of the c — M 
relation presented in this work. In addition, AppendixIClpro- 
vides a description of a simple code (available for download 
in IDL and Python) that computes concentrations, MAHs 
and accretion rates as a function of redshift for any cosmol¬ 
ogy- 


The c — M model presented in this work predicts a con¬ 
centration of c = 3 for microhalos of 10 - ' Mg at 2 = 31, 
in agreement with simula ti ons o f micr ohalo formation from 
lAnderhalden fe Diemand d2013l ) and fahiram 1 HoTl). In 
their work. lAnderhalden fc Diem andl ( 2013T) co mpare d the 
empirical c — M re lations from Bullock et al.l (1200 il l and 
iMaccio et al.l (l2008l l to their simulation outputs, and con¬ 
cluded that extrapolating simple power-law approximations 
to typical microhalo scales results in an overestimation of 
c by up to a factor of 10 at low—a. We also find large dif¬ 



Log 10 M [M s ] 


Figure 9. Dashed curves show the concentration—mass relation 
at various redshifts 2 . The solid and dashed lines show the c — M 
evolution of halos whose c(z) and M ( 2 ) are shown in the top and 
bottom panels of Fig. [8] respectively. All results were obtained 
using the physical model described in Sections [3] and [T] 


ferences between extrapolations of the high-mass power-law 
fits to low-masses and the predictions of our physical model, 
as can be seen by comparing the dashed and solid lines in 
Fig. [3 This impacts calculations of the dark matter (DM) 
annih ilation signal boost (see e.g. ISanchez-Conde fc Prada 
2014) and the power from DM annihilation (see e.g. [Mack 
2014T) . which make extensive use of c—M relations at various 
redshifts. In Section [ 6 ] we analyse this issue in more detail 
by calcula ting the pow er fro m DM annihilation, assuming 
either the I Duffy et al.l ( 2008 1 c—M relation (an empirical 
power-law fit for high halo masses) or the c — M relation 
from this work. 


5 EVOLUTION OF THE CONCENTRATION 

In this section we use our semi-analytic model to investigate 
the evolution of concentration and the effects that determine 
the slope of the c — M relation. The top panel of Fig.[ 8 ]shows 
the evolution of the concentration of halos that have masses 
of M 0 = 10 6 ,10 s , 10 10 ,10 12 and 10 14 M 0 at 2 = 0. The bot¬ 
tom panel shows the corresponding halo MAHs normalized 
to the final halo mass at 2 = 0 (Mo). We computed c(z) and 
M(z) following the models described in Sections l4l and 13721 
respectively. 

In Paper I, we used EPS theory to show that the 
MAH of all halos can be described by the expression 
M(z ) = Mo(l + z) a e ^ z , where the exponential is due to the 
fast growth at high —2 and the power-law due to the slow 
growth at low— 2 . In addition, the parameters a and j3 de¬ 
pend on halo mass. As a result, MAHs of halos larger than 
10 11 M 0 are mainly characterized by an exponential growth, 
whereas lower-mass halos exhibit a MAH closer to a power 
law, as can be seen in the bottom panel of Fig. [ 8 ] 

Comparing the coloured curves in the top and bot¬ 
tom panels of Fig. [ 8 ] we see an interesting relation between 
the evolution of c and the corresponding MAH. Dark mat¬ 
ter halos with a small growth rate are appear to contract, 
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and so their concentrations grow rapidly. This can be un¬ 
derstood as follows. At low redshift, during the dark en¬ 
ergy dominated epoch, M(z) of low-mass halos is charac¬ 
terized by a power law dCorrea et al.1 [2015al l. During this 
epoch, there is a drop in the accretion and merger rates of 
small halos, and the halo mass increases due to the evo¬ 
lution of the reference density used in the spherical over¬ 
density definition of the halo (p C iit(z) in this case). This 
so-called pseudo-evolution of the halo mass is thus driven 
by the halo mas s definition rather than the accretion of 
new material fsee lDiemer et af]|2013l and references therein). 
The pseudo-evolution of the halo mass gives the impression 
that concentratio ns are increasing because of contraction of 


the bound cor es (|Taai^ionii__el/aL 2004 j_ Zhao_e£_aL 2003|; 


iLu et al.l 120061 : iLi et al.l l2008i : Ivan den Bosch et al.l 20l4) . 
when in fact the core radius remain constant. Indeed, the 
evolution of the c—M relation has b een shown to be se nsitive 
to the definition of halo mass (e.g. iDuffv et al.ll200§l . If we 
assume that r_2 is constant in the redshift range z = 0 — 1, 
then the increase in r2oo due to the drop in p cr it gives the 
approximate increase in the concentration values. We find 

. 2 I ~ 1.4, in agreement 


c(z=0) _ r 2 oo(z=0) _ 

c(z = l) >’ 2 oo( z = 1 ) 


Pcrit (, _ 

, Pcrit ( z =0) J 


with the increase in concentration of a Mo = 10 6 M 0 halo. 

In the high-redshift regime (z !§> 1, matter-dominated 
epoch), the halo MAH is mainly characterized by exponen¬ 
tial growth. During this time, concentrations grow by a fac¬ 
tor of 2 (from z = 8 to z = 2) for a Mo = 10® Mq halo, 
decreasing to a factor of 1.08 (from 2 : = 8 to 2 = 2 ) for a 
Mo = 10 14 Mq halo. The pseudo-evolution of the halo mass 
is negligible in comparison with the high accretion rates, 
and the core radius increases simultaneously with the virial 
radius, hence the concentration hardly grows. 

In the case where the halo mass history is characterized 
by exponential growth at all z, representing the situation of 
a universe with no dark energy but fl m ^ 1 , we find that 
concentrations do not reach such large values at 3 = 0. We 
thus conclude that the evolution of the concentration is indi¬ 
rectly affected by the accelerated expansion of the Universe 
through the MAH and the halo mass definition. 

Next, we analyse how the evolution of the concentration 
determines the change in slope of the c—M relation. Fig. [9] 
shows c—M relations at various redshifts (dashed lines), and 
the c — M evolution of halos with Mo = 10®, 10 s , 10 10 ,10 12 
and 10 14 Mq (solid lines). From this figure, we see that the 
‘break’ in the low-redshift c — M relation that occurs at 
M ~ 10 11 Mq is produced by the change in the halo MAH. 

As mentioned, M(z) changes from being dominated by 
exponential growth for high-mass halos, to power-law growth 
for low-mass halos. It is natural to ask why the break in 
M(z) (and consequently in c — M) occurs at ~ 10 11 Mq. 
The answer is given by the rms of the linear theory den¬ 
sity perturbation field, a, that determines at which halo 
mass M(z) changes from power-law dominated to exponen¬ 
tial dominated. Since a > 1 at low masses (<§C 10 11 Mq) 
and a <C 1 at high masses (S> 10 11 Mq), the low values of a 
at high masses increase the value of the /3 parameter in the 
exponential function of the M(z) model, and thus M(z) is 
mostly dominated by the exponential growth. As a result, 
higher mass halos increase their mass faster, their inner cores 
increase with the virial radius, and their concentrations do 
not grow as rapidly. The different growth rate of the concen¬ 


trations produced by the change in the halos MAH, creates 
the ‘break’ in the c — M relation. 

Therefore, the break can be understood as being pro¬ 
duced by the varying power in the density perturbations 
through the halo MAH, where the MAH of low-mass halos 
at z < 1 is mostly driven by pseudo-evolution. The break is 
less prominent at z > 1, because at higher redshifts M(z) is 
mostly exponential for all halo masses (see eq. OH a —t 0 for 
Zi > 0 due to the growth factor and Fig. 0. 


6 IMPLICATIONS FOR THE DARK MATTER 
ANNIHILATION SIGNAL 

Dark matter (DM) particles are predicted to self-annihilate 
into Standard Model particl es, thus injecting ene rgy into 
the surrounding medium fe.g. lFurlanetto et al.ll200a l. In this 
section we calculate the DM annihilation rate per unit vol¬ 
ume produced by a smooth density field of DM (dominant 
before structure formation) and by cosmic structures (halos 
and microhalos). We separate DM into a smooth and struc¬ 
ture component because the spatial distribution of mass is 
almost completely smooth at very early cosmic times. Later 
gravitational instability causes ov erdensities to grow, until 
micro DM halos form. We follow ICirelli et al.l (12009(1 and 
[m ac 3 (1201414 in this calculation, and obtain the DM energy 
density and mean power from DM annihilation, assuming 
the physical c—M relation from this wo rk and from extrap¬ 
olations of the fits to simulations from IDuffv et all (120081) . 
Below we briefly describe the calculation of the DM annihi¬ 
lation rate produced by cosmic structures. 

The DM annihilation rate per unit volume results from 
the sum of two parts, a structure contribution and a smooth 
contribution. The smooth contribution, dominant before 
structure formation, z > 100 , can be written as 

R (A = t^j/Odm,o (1 + 2 ) ■ 

Where m x is the mass of the DM particle, (ou) the self- 
annihilation cross-section, and />dm,o the smooth DM den¬ 
sity today, pDM.O = HoMPcrit- 

The DM annihilation rate per unit volume due to halos, 
R, is given by 

= ^ J dr4nr 2 p 2 (r, M).(21) 

Here m x is the mass of the dark matter particle and ( av ) is 
the average annihilation cross section, which we assume to be 
100 G eV and 10 _ 2 ®cm 3 s _1 , respectively (e.g. lAprile fc et ahl 
1201214 . For the h alo mass function, j(z,M), we adopt the 
expression from lReed et al.l (120071 1. For p(r, M) we use the 
NFW density profile. We use M 200 as the halo mass defini¬ 
tion. 

Next, we calculate the effective DM density from struc¬ 
ture formation, defined as 


Pdm(z) 

= Pdm,o(1 + z) s Ri(z), 

(22) 


/2m 2 \ 1/2 


Ri(z) 

= (< 4 * w ) ■ 

(23) 

where 

Pdm,o is the average 

DM density today, 
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Figure 10. Left panel: Concentration-mass relation at various z (z = 0—10, top to bottom) from this work (blue solid lines) and from 
Duffy et al. (2008) (red dashed lines). At z = 10 and at a mass-scale of 1 Mq concentrations of Duffy et al. (2008) are a factor of 10 larger 
than concentrations predicted by this work, and a factor of 40 for a 10 —9 Mq halo. Middle panel: DM annihilation power per hydrogen 
nucleus as a function of redshift. The dashed line corresponds to the smooth component of the power and the dot-dashed lines to the 
structure co mponent. The solid lines show the sum of the two components (structure+smooth). The red lines correspond to the power 
assuming the iDuffv et al.1 (|2QQ8h c — M relation, whereas the blue lines correspond to the power assuming the c — M relation from this 
work. Right panel: as the middle panel, but showing the effective DM energy density as a function of redshift. 


£>dm,o = ^DMPcrit,o, with = 0.11. In this cal¬ 

culation we assume the WMAP5 cosmology in order to 
facilitate a comparison with models using the extrapolation 
of the Duffy et al. (2008) power-law fit to the c — M relation 
predicted by WMAP5 iV-body simulations. 

Similarly, we calculate the averaged volume power, per 
hydrogen nucleus, produced from DM annihilation events as 


P(z) = 2 m x c 2 LLL, 

n H (2j 

with n H (z) = n b p C rit,o(l 
density of hydrogen. 


(24) 

Yp)(l + 2) 3 /mH0 the number 


6.1 Implications 

Several of the models that h ave been used to predict the DM 
annihilation signal (see e .g. lLavalle et alj 120081 ; IPieri et al.l 
120081 : IPinzke et al.l 1201 ill have extrapolated c — M rela¬ 
tions, obtained from power-law fits to simulation results, 
to mass far below the resolution limit of the simulations. 
These power-law extrapolations assign huge concentrations 
to the smallest halos, thus increasing the DM annihilation 
power. In this section we explore how our physically moti¬ 
vated c — M relation, which flattens towards low-masses at 
low-redshift, affects the DM annihilation power when com- 
par ing it to the powe r calculated using an extrapolation of 
the lDuffv et al] (1200811 power-law c — M relation (the results 
are similar for other published power-law fits to the c — M 
relation from simulations). 

The left panel of Fig. [TO] shows a comparison between 
the concentration-mass relations at various re dshifts (z = 
0 — 10) from this work (blue solid lines) and of IDuffv et al.l 


5 f2 b is the present day baryon density parameter, Y p = 0.24 the 
primordial mass fraction of helium and raj the proton mass. 


d2008l l (red dashed lines). For z — 0 there is good agreement 
at the high-mass end between both relations. However, at 
2 — 1 0 and at a m ass-scale of lMg, the concentrations of 
IDuffv et al.l J2008ll are a factor of 10 larger than the concen¬ 
trations predicted by this work. For mass-scales of 10 -9 Mq 
the difference is a factor of 40. In the middle and right pan¬ 
els of Fig. [TO] we investigate the effects of these different 
c — M relations, showing the DM annihilation power per 
hydrogen nucleus and the effective DM energy density, re¬ 
spectively, as a function of redshift. In the middle panel, 
the dashed line corresponds to the smooth DM component 
of the power whereas the dot-dashed lines correspond to 
the structure component. The solid line shows the sum of 
the two components (structure+smooth). In e ach case the 
red lin es correspond to the power assuming the lDuffv et al] 
d2008li c — M relation, whereas the blue lines correspond to 
the power assuming the c — M relation from this work. The 
change in c—M affects the normalization of the power as 
well as the redshift at which structures begin to dominate. 
Lower concentrations result in lower central densities. Since 
the annihilation rate per unit volume, R, scales as p 2 , it is 
clear that R should decrease accordingly. At 2 = 0, the DM 
annihilation power that assumes the c—M relation predicted 
by this work is two orders of magnitude lower tha n the power 
obtained by extrapolating the lDuffv et al] (l2008h c—M rela¬ 
tion^ Ii+addition, the higher concentrations predicted by the 
IDuffv et al.l d200Sl f c — M relation imply that halos dominate 
the power over the smooth DM density component at higher 
redshifts. Adopting the c—M relation from this work results 
in the power from structures starting to do minate at 2 m 50 
rather than at 2 » 85 (in agreement with iMac 3 120141 and 
iNg et af1l2014lh This lower redshift of structure formation 
dominating over the smooth component could have impor¬ 
tant implicatio ns for searches of the ‘D ark Ages’ by radio 
telescopes (e.g. IPritchard fe Loelj|20T3 ). 
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Figure 11. Distribution of particles within r 200 - We show the 
particle distribution of two different halo samples, the first sample 
contains halos of 10 15 Mq identified at z% = 0 (solid lines) and 
the second sample contains halos of 10 14 Mq identified at Zi = 1 
(dashed lines). The solid grey line corresponds to a NFW density 
profile with concentration c = 4 and the vertical dashed line to 
the corresponding r _2 radius. The different colours indicate the 
redshift (z a ccr) during which the particles were accreted on to the 
halos (note that the same colour corresponds to different redshift 
ranges for different Zj). 


7 DISCUSSION 

Dark matter halo concentrations have recently been the 
subject of extensive analysis due to the controversial re¬ 
ports of an upturn at the high-mass end of the relation 
(1 Munoz-Cuart as et al.ll201ll:IPrada et al]l2012l : lKlvpin et alJ 

l2014l ; lDiemer &; Kravtsov! 20 id ). However, the semi-analytic 

model for dark matter halo concentrations presented in this 
work does not predict such an upturn. In this section, we 
review the main assumptions that the model relies on and 
discuss the plausibility of the existence of the upturn. 

First, the model assumes that the halo density profile 
is described by the NFW profile at all ti mes. A lthough it is 
known that the Einasto profile dEinastol Il965h is more ac- 
curate than the NFW profile dGao et al.ll200a ). it has an 
extra ‘shape’ parameter that complicates the fitting proce¬ 
dure and affects the concentration. Also, the residuals from 
the systematic deviations from the NFW shape are gener¬ 
ally smaller than 10% and the NFW concentrations only 
differ by 10 — 20% from Einas to fit and the velocity pro¬ 
file fit (as recently discussed bv lDutton fc Maccioll2014l and 
iKlypin et al.l 12014 1. We thus conclude that using the NFW 
profile to predict densities is not a major determinant in the 
model. 

Secondly, the model depends on the calibration of the 
Pcrit( 2 —2) — ( p)(< r~ 2) relation, which implies that halo for¬ 
mation is an ‘inside out’ process, where the central part 
of a dark matter halo (contained within r- 2 ) forms Erst, 
and later accretion and mergers increase the mass and size 
of the halo without adding much material to its inner re¬ 
gions tHuss et al ] 19991 : IWang fc White! 120091 9 . We test the 
assumption of ‘halo formation’ of our model. We analyse 
the distribution of particles within r 2 oo, and differentiate the 


particles according to the period of time during which they 
were accreted. Fig. Illl shows the radial distribution of parti¬ 
cles of two different halo samples. The first sample contains 
halos of 10 15 Mq identified at Zi = 0 (which are formed at 
Z -2 ~ 1) and the second sample contains halos of 10 14 Mq 
identified at Zi = 1 (formed at Z -2 ~ 2). We analysed halos 
of different masses but focused on massive halos, because 
their large radial velocities makes them more likely to con¬ 
tain recently accreted particles in their inner regions, and 
they are located in the upturn in the c — M relation. 


Fig. [TT] shows the NFW density profile (in grey solid 
line) and the r _2 radius (in vertical dashed line) for a con¬ 
centration of c = 4. The figure also shows the distribution 
of particles at Zi = 0 (solid coloured lines) and at Z; = 1 
(dashed coloured lines). The different colours indicate the 
redshift (z acc r) during which the particles were accreted on 
to the halos. The blue lines show that recently accreted par¬ 
ticles are distributed around T 2 oo and that only a tiny frac¬ 
tion (< 2 % of the total, in the two cases), reside in the inner 
parts of the halo. In the case of the distribution of particles 
at Zi = 0, 2.5% of the total particles are in regions within 
r -2 after being accreted during z acc r = 0.25 — 0.5, and 4.5% 
during z acC r = 0.5 — 1. The same behaviour is observed in ha¬ 
los of different masses identified at higher redshifts. We find 
that 8 % of all particles accreted after the halo has formed 
are in the centre, not enough to significantly alter the mass 
within r -2 so as to increase concentration. We then find the 
halo formation assumption that the model relies on to be 
valid. 


Finally, in the calibration of the p C rit(z- 2 ) — (p)(< r- 2 ) 
relation, we only consider relaxed halos. The selection con¬ 
ditions g enerally used to dif f erentiate relaxed halos from un¬ 
relaxed (iMaccio et al.l 20071: Neto et al.ll 2007h have recently 
been revisited by Klypin et al.l ( 20141 '). These conditions in¬ 
clude the virial parameters [2K/\W\ — 1, where K and W are 
the kinetic and potential energies), the offset parameter X 0 h 
(distance between the potential minimum and the center of 
mass) , and the spin parameter. In their work, iKlypin et al.l 
(l2014h argued that the virial equilibrium condition is too 
simplistic and needs to include the effects of the surface 
pressure and external forces. They applied these corrections 
to the virial parameters and selected halos that had previ¬ 
ously been rejected. As a result, they obtained an upturn in 
the high-mass end of the c — M relation and claimed that 
the large concentration of massive halos is due to their in¬ 
fall velocities, which are more radial and result in deeper 
penetration of infalling mass into the halo that reaches the 
inner parts. In this work we selected relaxed h alos usin g only 
the condition that X 0 a < 0.07 following iDuffv et al.l (120080 
and lNeto et al.l (12007 1. who found that this simple criterion 
resulted in the removal of the vast majority of unrelaxed 
haloes. We did not use any additional criteria and did not 
find any upturn at high halo masses, but concluded that the 
strong flattening of the c — M relation at high redshift is due 
to unrelaxed h alos (Fig. El). We cannot say we disagree with 
IKlypin et al.l (1201419 regarding the relaxation conditions, be¬ 
cause our simulations do not have sufficiently large box sizes 
to model a large sample of the massive highly unrelaxed ha¬ 
los that likely ‘shape’ the upturn. 
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8 CONCLUSION 

In this paper, we have linked the concentration of a halo to 
its MAH. We extended the analytic framework presented in 
Paper I to show that the halo mass history M(z, M{z/), z/) 
of a halo with mass M{z/) at Zi can be described by 


M(z,M(zi), Zi) = M(zi)( 1 + z — Zi) a e' 




where a and /3 are parameters that depend on M(zi) and 
cosmology. We have compared the above formula to simu¬ 
lation outputs and the mo st recent empirical mass hi story 
model from the literature (Ivan den Bosch et al] 12014! ) and 
found generally reasonable ag reement. _ _ 

Building on the work bv iLudlow et al.1 120141 ). we pre¬ 
sented a semi-analytic model for halo concentration as a 
function of halo mass and redshift. The model connects the 
analytic model for the MAH from Paper I to halo concen¬ 
tration through an empirical relation between concentra¬ 
tion and formation redshift, obtained through fits to simu- 
lation data in Paper II. The formation redsh ift definition of 
ILudlow et al.l J2013f) and lLudlow et al.ll20l4 2 _ 2 , defined as 
the redshift at which the mass of the main progenitor equals 
the mass enclosed within the scale radius of the NFW den¬ 
sity profile at z = Zi, results in an inner halo structure that 
reflects the background density of the Universe at the time 
when the halo formed. 

The resulting c — M relations were tested using N- 
body simulations, and compared to th e most recent empiri¬ 
cal c — M relations from t he lit e rature (Ivan den Bosch et al.l 
120141 : iDutton fe Maccidl 120141 ; IDiemer fe Kravtsov! 120151) . 
The ‘upturn ’ at h igh masses seen by some studies 
JlGvmn_et_al. 20 111: Prada et all 120121 : IDutton fe Maccidl 
120141 : IDiemer fe Kravtsovl 1 2(Jl sT T is not reproduced by our 
physically derived model which, however, only applies to re¬ 
laxed halos. We analysed the c—M relations obtained from 
the numerical simulations, where we differentiated between 
relaxed and unrelaxed halos. We found that the upturn is 
due to the inc lusion of unrelaxed halos, supporting the pre¬ 
vious claim of lLudlow et al.l (120121) . 

We applied our model to a large range of halo mass 
(log 10 M/Mq = [—2,16]) and redshift (2 = 0 — 20), and 
provided fits to the c — M relations as well as numerical 
routine^ to compute concentrations and MAHs as a func¬ 
tion of halo mass, redshift and cosmology. We caution the 
reader that baryonic processes will almost certainly increase 
the scatter in the c — M relation and will modify the inner 
density profile (e.g.lDuffv et al.ll2 010l: [Govcrnato et a l 1120121: 


iTevssier et al.ll2013: Schaller et al. 20151') . 

Our model predicts a change in the slope of the c — M 
relation at z = 0 — 3 and a ‘break’ in the z ~ 0 c—M relation 
at a mass of ~ 10 11 Mg. We analysed the evolution of con¬ 
centration and found that it increases more rapidly during 
the dark energy era, when the accretion rates of dark matter 
halos decrease due to the accelerated expansion of the Uni¬ 
verse. We found that the break at a halo mass ~ 10 11 M@ 
results from the change of the functional form of M(z), that 


6 Available at https://bitbucket.org/astroduff/commah and 
http://astro.physics.unimelb.edu.au/ in Research/Public-Data- 
Releases/COMMAH. See Appendix ICl for a short description. 


goes from being dominated by a power-law (for low-mass ha¬ 
los) to an exponential (for high-mass halos). This change in 
M(z) is driven by the rms of the linear density perturbation 
field. Halos with mass M 10 11 Mg, have cr 1, are char¬ 
acterized by an exponential growth, and have lower concen¬ 
trations as a result. Whereas halos with mass M 10 11 M©, 
have large < 7 , are characterized by a power-law growth. In 
this last case, there is a pseudo-evolution in the halo masses 
(i.e. mass growth due to t he definition of the halo in terms 
of an overdensity criterion, IDiemer et al . IB) and the core 
radius remains approximately constant, causing the concen¬ 
trations to grow. The different growth rate of the concentra¬ 
tions at low and high mass produces the break in the 2^0 
c—M relation. This break is not so evident for z > 1, because 
at higher redshifts M(z) is mostly exponential for all halo 
masses (a —> 0 for Zi > 0 due to the growth factor), causing 
all concentrations to grow at approximately the same rate 
(as seen in Fig. El). 

Finally, we addressed the impact of the c — M relation 
presented in this work on predictions for the dark matter 
annihilation signal. We calculated the DM annihilation rate 
from cosmic struct ures and compared the results obtained 
by extrapolating the lDuffv et all (120081 ) power-law fit (which 
is similar to other published fits to the results of simulations) 
to the rate obtained by using the c — M relation predicted by 
our model. We found that the power from DM annihilation 
at 2 = 0 is two orders of magnitude lower than the power 
obtained by extrapolating the Duffy et al. (2008) c — M 
relation (in agreement with IUackll2014l : iNg et al 1 12014) . 
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APPENDIX A: SIMULATIONS AND 
COSMOLOGY 

In this work, we use the set of cos mological DMONLY sim¬ 
ulations from the OWLS project (ISchave et alJlioiol h The 
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Table Al. List of simulations. From left-to-right the columns show: simulation identifier; comoving box size; number of dark matter 
particles; dark matter particle mass; comoving (Plummer-equivalent) gravitational softening; maximum physical softening. 


Simulation 

L 

(h — 1 Mpc) 

N 

m dm 

r'M 0 ) 

^com 

(/i _1 kpc) 

£prop 

(h~ ^Tpc) 

DMONLY_WMAP5_L400N512 

400 

512 3 

3.4 x 10 10 

31.25 

8.00 

DMONLY_WMAP5_L200N512 

200 

512 3 

3.2 x 10 9 

15.62 

4.00 

DMONLY_WMAP5_L100N512 

100 

512 3 

5.3 x 10 8 

7.81 

2.00 

DMONLY_WMAP5_L050N512 

50 

512 3 

6.1 x 10 7 

3.91 

1.00 

DMONLY_WMAP5_L025N512 

25 

512 3 

8.3 x 10 6 

2.00 

0.50 


Table A2. Cosmological parameters. 


Simulation 


Pa 

h 

^8 

n s 

DMONLY.WMAPl 

0.25 

0.75 

0.73 

0.90 

1.000 

DMONLY-WMAP3 

0.238 

0.762 

0.73 

0.74 

0.951 

DMONLY_WMAP5 

0.258 

0.742 

0.72 

0.796 

0.963 

DMONLY-WMAP9 

0.282 

0.718 

0.70 

0.817 

0.964 

DMONLY-Planckl 

0.317 

0.683 

0.67 

0.834 

0.962 


particle masses and gravitational softenings for each of these 
simulations are listed in Table I All Table IA2I lists the sets 
of cosmological parameters adopted in the different simula¬ 
tions. In this section we describe the numerical techniques 
used to compute merger trees. 

The first step towards studying the mass assembly his¬ 
tory of halos is to identify gravitationally bound structures 
and build halo merger trees. We begin by selec ting the 
largest halo in each FoF group dDavis et al.l 1 19851 : we use 
a linking length of 0.2) (i.e. the main subhalo of FoF groups 
that is not embedded inside larger halos). Halo virial masses 
and radii are determined using a sphe rical overdensity rou - 
tine within the SUBFIND algorithm (ISpringel et aklhiOOlh 
centred on the main subhalo of FoF halos. For a complete 
description on the method used to build the halo merger 
trees, see Paper II. 


APPENDIX B: COSMOLOGY DEPENDENCE 

The adopted cosmological parameters affect the halo MAH 
so that the larger ag or n m , the earlier halos assemble. As the 
formation time increases with increasing ag or so does 
concentration. Therefore, in this section, we analyse how the 
change in cosmology affects our c — M relation model. 

Our c—M model described in Section[4]relies on the halo 
mass accretion history model, which we showed in Section [I] 
to be suitable for any cosmology. However, our semi-analytic 
model for halo concentration also relies on the formation 
redshift through the best-fitting relation given by eq. 113, 

(p)(< r~ 2 ,Zj) = A x Pcrit(g —2) 

/2(:ri 1 ('S-,;) Peril (Ur ) 

where A depends on cosmology (A = 887 for WMAP5). 

We investigate the cosmology dependence of A by fol¬ 
lowing the analysis done in Section 4.1 and using the simu¬ 
lations with different cosmologies listed in Table 1X21 We cal¬ 
culate the best-fitting (p)(< r_ 2 ,z;) — p C rit(«- 2 , zl) relation 


to obtain the parameter A cosmo , where cosmo is WMAP1, 
WMAP3, WMAP9 or Planck. We found that if we keep 
A = 887 fixed when computing the concentration-mass re¬ 
lations for cosmologies other than WMAP5, we obtain rela¬ 
tions that are in very good agreement with the various rela¬ 
tions in the literature. However, we obtain even better agree¬ 
ment with the simulations when we let A vary slightly with 
cosmology. We found that Awmapi = 853, Awmap3 = 850, 
Awmap9 = 950 and A P i anc k = 880, reproduce the c—M 
relations best. For a calculation of c—M relations in cos¬ 
mologies other than the ones listed above, we recommend 
the reader to set A = 887 fixed. 

It is important to note that if the A parameter is 
in fact cosmology dependent then halos which formed at 
different formation times (e.g. z_2,cosmoi and z_2,cosmo2), 
but that correspond to the same background density 
(Pcrit (^— 2 ,cosmol ) = Pcrit (-Z—2,cosmo2)), will have different con¬ 
centrations (Ccosmoi 7 ^ c C osmo 2 ) ■ This implies that other fac¬ 
tors affect the dark matter halo density profiles, and that 
the halo MAH alone is not sufficient for predicting concen¬ 
trations. 


B1 Fitting functions for the c — M relation 

In this section, we provide fitting functions for the c—M 
relation in the high-z and low-z regimes for Planck cosmol¬ 
ogy. The following expression is suitable for the low-redshift 
regime (z ^ 4) and at all halo masses, 


logi 0 c = a + P log 10 (M/ M 0 )[l + 7 (log 10 M/ M 0 ) 2 ], 
a = 1.7543 - 0.2766(1 + z) + 0.02039(1 +z) 2 , 

P = 0.2753 + 0.00351(1 + z) - 0.3038(1 +z) 0 0269 , 
7 = -0.01537 + 0.02102(1 +z)~ 0 1475 . 


In the high-redshift regime the c—M relation can be 
fitted using only two parameters. The following expression 
is suitable for z > 4 and at all halo masses, 


l°gio c = ex + P log 10 (M/ M 0 ), 

a = 1.3081 - 0.1078(1 + z) +0.00398(1 +z) 2 , 
P = 0.0223 - 0.0944(1 + z) -0 ' 3907 . 


APPENDIX C: ONLINE MATERIAL 

The concentration-mass relation model presented in this 
work, as well as the halo mass accretion history 
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model, are included in the code named COMMAH for 
COncentration-Mass relation and Mass Accretion History, 
available at https://bitbucket.org/astroduff/commah and 
http://astro.physics.unimelb.edu.au/ in Research/Public- 
Data-Releases/COMMAH. The code is available in both 
the python and IDL languages. Also, COMMAH is in the PyPi 
python package, to install it type ‘pip install commah’. In 
this section we present a short overview of COMMAH. 

COMMAH is a routine that follows the analytic model de¬ 
scribed in Section [3] to calculate the MAH of a halo of mass 
Mo at z = 0 in any given redshift interval (e.g. M(z) between 
2 = 0 and 10). Also, COMMAH calculates halo concentrations 
following the semi-analytic model described in Section[f] and 
outputs the c — M relation at any given redshift. In addition, 
it also computes the dark matter accretion rate, the rms of 
the density field, peak height, and the integral of the NFW 
density profile (see eq. ED, suitable for DM annihilation cal¬ 
culations. COMMAH calculates concentration solving eqns. m 
and m by performing a Levenberg-Marquardt method. As 
described in the previous section, commah is suitable for any 
cosmology. 


